3.3.40 \(\int \frac {x^2}{(d+e x^2) (a+c x^4)} \, dx\) [240]
Optimal. Leaf size=337 \[ -\frac {\sqrt {d} \sqrt {e} \tan ^{-1}\left (\frac {\sqrt {e} x}{\sqrt {d}}\right )}{c d^2+a e^2}-\frac {\left (\sqrt {c} d+\sqrt {a} e\right ) \tan ^{-1}\left (1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt {2} \sqrt [4]{a} \sqrt [4]{c} \left (c d^2+a e^2\right )}+\frac {\left (\sqrt {c} d+\sqrt {a} e\right ) \tan ^{-1}\left (1+\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt {2} \sqrt [4]{a} \sqrt [4]{c} \left (c d^2+a e^2\right )}+\frac {\left (\sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} \sqrt [4]{c} \left (c d^2+a e^2\right )}-\frac {\left (\sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} \sqrt [4]{c} \left (c d^2+a e^2\right )} \]
[Out]
1/8*ln(-a^(1/4)*c^(1/4)*x*2^(1/2)+a^(1/2)+x^2*c^(1/2))*(-e*a^(1/2)+d*c^(1/2))/a^(1/4)/c^(1/4)/(a*e^2+c*d^2)*2^
(1/2)-1/8*ln(a^(1/4)*c^(1/4)*x*2^(1/2)+a^(1/2)+x^2*c^(1/2))*(-e*a^(1/2)+d*c^(1/2))/a^(1/4)/c^(1/4)/(a*e^2+c*d^
2)*2^(1/2)+1/4*arctan(-1+c^(1/4)*x*2^(1/2)/a^(1/4))*(e*a^(1/2)+d*c^(1/2))/a^(1/4)/c^(1/4)/(a*e^2+c*d^2)*2^(1/2
)+1/4*arctan(1+c^(1/4)*x*2^(1/2)/a^(1/4))*(e*a^(1/2)+d*c^(1/2))/a^(1/4)/c^(1/4)/(a*e^2+c*d^2)*2^(1/2)-arctan(x
*e^(1/2)/d^(1/2))*d^(1/2)*e^(1/2)/(a*e^2+c*d^2)
________________________________________________________________________________________
Rubi [A]
time = 0.17, antiderivative size = 337, normalized size of antiderivative = 1.00, number of steps
used = 12, number of rules used = 8, integrand size = 22, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.364, Rules used = {1302, 211,
1182, 1176, 631, 210, 1179, 642} \begin {gather*} -\frac {\sqrt {d} \sqrt {e} \text {ArcTan}\left (\frac {\sqrt {e} x}{\sqrt {d}}\right )}{a e^2+c d^2}-\frac {\text {ArcTan}\left (1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right ) \left (\sqrt {a} e+\sqrt {c} d\right )}{2 \sqrt {2} \sqrt [4]{a} \sqrt [4]{c} \left (a e^2+c d^2\right )}+\frac {\text {ArcTan}\left (\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}+1\right ) \left (\sqrt {a} e+\sqrt {c} d\right )}{2 \sqrt {2} \sqrt [4]{a} \sqrt [4]{c} \left (a e^2+c d^2\right )}+\frac {\left (\sqrt {c} d-\sqrt {a} e\right ) \log \left (-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {a}+\sqrt {c} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} \sqrt [4]{c} \left (a e^2+c d^2\right )}-\frac {\left (\sqrt {c} d-\sqrt {a} e\right ) \log \left (\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {a}+\sqrt {c} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} \sqrt [4]{c} \left (a e^2+c d^2\right )} \end {gather*}
Antiderivative was successfully verified.
[In]
Int[x^2/((d + e*x^2)*(a + c*x^4)),x]
[Out]
-((Sqrt[d]*Sqrt[e]*ArcTan[(Sqrt[e]*x)/Sqrt[d]])/(c*d^2 + a*e^2)) - ((Sqrt[c]*d + Sqrt[a]*e)*ArcTan[1 - (Sqrt[2
]*c^(1/4)*x)/a^(1/4)])/(2*Sqrt[2]*a^(1/4)*c^(1/4)*(c*d^2 + a*e^2)) + ((Sqrt[c]*d + Sqrt[a]*e)*ArcTan[1 + (Sqrt
[2]*c^(1/4)*x)/a^(1/4)])/(2*Sqrt[2]*a^(1/4)*c^(1/4)*(c*d^2 + a*e^2)) + ((Sqrt[c]*d - Sqrt[a]*e)*Log[Sqrt[a] -
Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/(4*Sqrt[2]*a^(1/4)*c^(1/4)*(c*d^2 + a*e^2)) - ((Sqrt[c]*d - Sqrt[a]*
e)*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2])/(4*Sqrt[2]*a^(1/4)*c^(1/4)*(c*d^2 + a*e^2))
Rule 210
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(-(Rt[-a, 2]*Rt[-b, 2])^(-1))*ArcTan[Rt[-b, 2]*(x/Rt[-a, 2])
], x] /; FreeQ[{a, b}, x] && PosQ[a/b] && (LtQ[a, 0] || LtQ[b, 0])
Rule 211
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]
Rule 631
Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> With[{q = 1 - 4*Simplify[a*(c/b^2)]}, Dist[-2/b, Sub
st[Int[1/(q - x^2), x], x, 1 + 2*c*(x/b)], x] /; RationalQ[q] && (EqQ[q^2, 1] || !RationalQ[b^2 - 4*a*c])] /;
FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]
Rule 642
Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[d*(Log[RemoveContent[a + b*x +
c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]
Rule 1176
Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[2*(d/e), 2]}, Dist[e/(2*c), Int[1/S
imp[d/e + q*x + x^2, x], x], x] + Dist[e/(2*c), Int[1/Simp[d/e - q*x + x^2, x], x], x]] /; FreeQ[{a, c, d, e},
x] && EqQ[c*d^2 - a*e^2, 0] && PosQ[d*e]
Rule 1179
Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[-2*(d/e), 2]}, Dist[e/(2*c*q), Int[
(q - 2*x)/Simp[d/e + q*x - x^2, x], x], x] + Dist[e/(2*c*q), Int[(q + 2*x)/Simp[d/e - q*x - x^2, x], x], x]] /
; FreeQ[{a, c, d, e}, x] && EqQ[c*d^2 - a*e^2, 0] && NegQ[d*e]
Rule 1182
Int[((d_) + (e_.)*(x_)^2)/((a_) + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[a*c, 2]}, Dist[(d*q + a*e)/(2*a*c),
Int[(q + c*x^2)/(a + c*x^4), x], x] + Dist[(d*q - a*e)/(2*a*c), Int[(q - c*x^2)/(a + c*x^4), x], x]] /; FreeQ
[{a, c, d, e}, x] && NeQ[c*d^2 + a*e^2, 0] && NeQ[c*d^2 - a*e^2, 0] && NegQ[(-a)*c]
Rule 1302
Int[(((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.))/((a_) + (c_.)*(x_)^4), x_Symbol] :> Int[ExpandIntegrand[(
f*x)^m*((d + e*x^2)^q/(a + c*x^4)), x], x] /; FreeQ[{a, c, d, e, f, m}, x] && IntegerQ[q] && IntegerQ[m]
Rubi steps
\begin {align*} \int \frac {x^2}{\left (d+e x^2\right ) \left (a+c x^4\right )} \, dx &=\int \left (-\frac {d e}{\left (c d^2+a e^2\right ) \left (d+e x^2\right )}+\frac {a e+c d x^2}{\left (c d^2+a e^2\right ) \left (a+c x^4\right )}\right ) \, dx\\ &=\frac {\int \frac {a e+c d x^2}{a+c x^4} \, dx}{c d^2+a e^2}-\frac {(d e) \int \frac {1}{d+e x^2} \, dx}{c d^2+a e^2}\\ &=-\frac {\sqrt {d} \sqrt {e} \tan ^{-1}\left (\frac {\sqrt {e} x}{\sqrt {d}}\right )}{c d^2+a e^2}-\frac {\left (d-\frac {\sqrt {a} e}{\sqrt {c}}\right ) \int \frac {\sqrt {a} \sqrt {c}-c x^2}{a+c x^4} \, dx}{2 \left (c d^2+a e^2\right )}+\frac {\left (d+\frac {\sqrt {a} e}{\sqrt {c}}\right ) \int \frac {\sqrt {a} \sqrt {c}+c x^2}{a+c x^4} \, dx}{2 \left (c d^2+a e^2\right )}\\ &=-\frac {\sqrt {d} \sqrt {e} \tan ^{-1}\left (\frac {\sqrt {e} x}{\sqrt {d}}\right )}{c d^2+a e^2}+\frac {\left (\sqrt [4]{c} \left (d-\frac {\sqrt {a} e}{\sqrt {c}}\right )\right ) \int \frac {\frac {\sqrt {2} \sqrt [4]{a}}{\sqrt [4]{c}}+2 x}{-\frac {\sqrt {a}}{\sqrt {c}}-\frac {\sqrt {2} \sqrt [4]{a} x}{\sqrt [4]{c}}-x^2} \, dx}{4 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}+\frac {\left (\sqrt [4]{c} \left (d-\frac {\sqrt {a} e}{\sqrt {c}}\right )\right ) \int \frac {\frac {\sqrt {2} \sqrt [4]{a}}{\sqrt [4]{c}}-2 x}{-\frac {\sqrt {a}}{\sqrt {c}}+\frac {\sqrt {2} \sqrt [4]{a} x}{\sqrt [4]{c}}-x^2} \, dx}{4 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}+\frac {\left (d+\frac {\sqrt {a} e}{\sqrt {c}}\right ) \int \frac {1}{\frac {\sqrt {a}}{\sqrt {c}}-\frac {\sqrt {2} \sqrt [4]{a} x}{\sqrt [4]{c}}+x^2} \, dx}{4 \left (c d^2+a e^2\right )}+\frac {\left (d+\frac {\sqrt {a} e}{\sqrt {c}}\right ) \int \frac {1}{\frac {\sqrt {a}}{\sqrt {c}}+\frac {\sqrt {2} \sqrt [4]{a} x}{\sqrt [4]{c}}+x^2} \, dx}{4 \left (c d^2+a e^2\right )}\\ &=-\frac {\sqrt {d} \sqrt {e} \tan ^{-1}\left (\frac {\sqrt {e} x}{\sqrt {d}}\right )}{c d^2+a e^2}+\frac {\sqrt [4]{c} \left (d-\frac {\sqrt {a} e}{\sqrt {c}}\right ) \log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}-\frac {\sqrt [4]{c} \left (d-\frac {\sqrt {a} e}{\sqrt {c}}\right ) \log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}+\frac {\left (\sqrt [4]{c} \left (d+\frac {\sqrt {a} e}{\sqrt {c}}\right )\right ) \text {Subst}\left (\int \frac {1}{-1-x^2} \, dx,x,1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}-\frac {\left (\sqrt [4]{c} \left (d+\frac {\sqrt {a} e}{\sqrt {c}}\right )\right ) \text {Subst}\left (\int \frac {1}{-1-x^2} \, dx,x,1+\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}\\ &=-\frac {\sqrt {d} \sqrt {e} \tan ^{-1}\left (\frac {\sqrt {e} x}{\sqrt {d}}\right )}{c d^2+a e^2}-\frac {\sqrt [4]{c} \left (d+\frac {\sqrt {a} e}{\sqrt {c}}\right ) \tan ^{-1}\left (1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}+\frac {\sqrt [4]{c} \left (d+\frac {\sqrt {a} e}{\sqrt {c}}\right ) \tan ^{-1}\left (1+\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )}{2 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}+\frac {\sqrt [4]{c} \left (d-\frac {\sqrt {a} e}{\sqrt {c}}\right ) \log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}-\frac {\sqrt [4]{c} \left (d-\frac {\sqrt {a} e}{\sqrt {c}}\right ) \log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} \left (c d^2+a e^2\right )}\\ \end {align*}
________________________________________________________________________________________
Mathematica [A]
time = 0.08, size = 232, normalized size = 0.69 \begin {gather*} \frac {-8 \sqrt [4]{a} \sqrt [4]{c} \sqrt {d} \sqrt {e} \tan ^{-1}\left (\frac {\sqrt {e} x}{\sqrt {d}}\right )+\sqrt {2} \left (-2 \left (\sqrt {c} d+\sqrt {a} e\right ) \tan ^{-1}\left (1-\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )+2 \left (\sqrt {c} d+\sqrt {a} e\right ) \tan ^{-1}\left (1+\frac {\sqrt {2} \sqrt [4]{c} x}{\sqrt [4]{a}}\right )+\left (\sqrt {c} d-\sqrt {a} e\right ) \left (\log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )-\log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{c} x+\sqrt {c} x^2\right )\right )\right )}{8 \sqrt [4]{a} \sqrt [4]{c} \left (c d^2+a e^2\right )} \end {gather*}
Antiderivative was successfully verified.
[In]
Integrate[x^2/((d + e*x^2)*(a + c*x^4)),x]
[Out]
(-8*a^(1/4)*c^(1/4)*Sqrt[d]*Sqrt[e]*ArcTan[(Sqrt[e]*x)/Sqrt[d]] + Sqrt[2]*(-2*(Sqrt[c]*d + Sqrt[a]*e)*ArcTan[1
- (Sqrt[2]*c^(1/4)*x)/a^(1/4)] + 2*(Sqrt[c]*d + Sqrt[a]*e)*ArcTan[1 + (Sqrt[2]*c^(1/4)*x)/a^(1/4)] + (Sqrt[c]
*d - Sqrt[a]*e)*(Log[Sqrt[a] - Sqrt[2]*a^(1/4)*c^(1/4)*x + Sqrt[c]*x^2] - Log[Sqrt[a] + Sqrt[2]*a^(1/4)*c^(1/4
)*x + Sqrt[c]*x^2])))/(8*a^(1/4)*c^(1/4)*(c*d^2 + a*e^2))
________________________________________________________________________________________
Maple [A]
time = 0.18, size = 246, normalized size = 0.73
| | |
method |
result |
size |
| | |
default |
\(-\frac {d e \arctan \left (\frac {e x}{\sqrt {d e}}\right )}{\left (a \,e^{2}+c \,d^{2}\right ) \sqrt {d e}}+\frac {\frac {e \left (\frac {a}{c}\right )^{\frac {1}{4}} \sqrt {2}\, \left (\ln \left (\frac {x^{2}+\left (\frac {a}{c}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{c}}}{x^{2}-\left (\frac {a}{c}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{c}}}\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{c}\right )^{\frac {1}{4}}}+1\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{c}\right )^{\frac {1}{4}}}-1\right )\right )}{8}+\frac {d \sqrt {2}\, \left (\ln \left (\frac {x^{2}-\left (\frac {a}{c}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{c}}}{x^{2}+\left (\frac {a}{c}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{c}}}\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{c}\right )^{\frac {1}{4}}}+1\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{c}\right )^{\frac {1}{4}}}-1\right )\right )}{8 \left (\frac {a}{c}\right )^{\frac {1}{4}}}}{a \,e^{2}+c \,d^{2}}\) |
\(246\) |
risch |
\(\frac {\sqrt {-d e}\, \ln \left (\left (-16 \left (-d e \right )^{\frac {5}{2}} a^{2} c \,e^{3}+16 \left (-d e \right )^{\frac {5}{2}} a \,c^{2} d^{2} e -14 a^{2} c d \,e^{4} \left (-d e \right )^{\frac {3}{2}}+20 a \,c^{2} d^{3} e^{2} \left (-d e \right )^{\frac {3}{2}}+2 d^{5} \left (-d e \right )^{\frac {3}{2}} c^{3}-\sqrt {-d e}\, a^{3} e^{7}+3 \sqrt {-d e}\, a \,c^{2} d^{4} e^{3}+2 \sqrt {-d e}\, c^{3} d^{6} e \right ) x -a^{3} d \,e^{7}-2 a^{2} c \,d^{3} e^{5}-a \,c^{2} d^{5} e^{3}\right )}{2 a \,e^{2}+2 c \,d^{2}}-\frac {\sqrt {-d e}\, \ln \left (\left (16 \left (-d e \right )^{\frac {5}{2}} a^{2} c \,e^{3}-16 \left (-d e \right )^{\frac {5}{2}} a \,c^{2} d^{2} e +14 a^{2} c d \,e^{4} \left (-d e \right )^{\frac {3}{2}}-20 a \,c^{2} d^{3} e^{2} \left (-d e \right )^{\frac {3}{2}}-2 d^{5} \left (-d e \right )^{\frac {3}{2}} c^{3}+\sqrt {-d e}\, a^{3} e^{7}-3 \sqrt {-d e}\, a \,c^{2} d^{4} e^{3}-2 \sqrt {-d e}\, c^{3} d^{6} e \right ) x -a^{3} d \,e^{7}-2 a^{2} c \,d^{3} e^{5}-a \,c^{2} d^{5} e^{3}\right )}{2 \left (a \,e^{2}+c \,d^{2}\right )}+\frac {\left (\munderset {\textit {\_R} =\RootOf \left (1+\left (a^{3} c \,e^{4}+2 a^{2} c^{2} d^{2} e^{2}+c^{3} a \,d^{4}\right ) \textit {\_Z}^{4}+4 a c d e \,\textit {\_Z}^{2}\right )}{\sum }\textit {\_R} \ln \left (\left (\left (-2 a^{4} c \,e^{7}-2 a^{3} c^{2} d^{2} e^{5}+2 a^{2} c^{3} d^{4} e^{3}+2 a \,c^{4} d^{6} e \right ) \textit {\_R}^{5}+\left (-7 a^{2} c d \,e^{4}+10 a \,c^{2} d^{3} e^{2}+c^{3} d^{5}\right ) \textit {\_R}^{3}+\left (-2 a \,e^{3}+4 c \,d^{2} e \right ) \textit {\_R} \right ) x +\left (-3 a^{3} c d \,e^{5}-6 a^{2} c^{2} d^{3} e^{3}-3 c^{3} d^{5} e a \right ) \textit {\_R}^{4}+\left (-13 a c \,d^{2} e^{2}-c^{2} d^{4}\right ) \textit {\_R}^{2}-4 d e \right )\right )}{4}\) |
\(598\) |
| | |
|
|
|
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(x^2/(e*x^2+d)/(c*x^4+a),x,method=_RETURNVERBOSE)
[Out]
-d*e/(a*e^2+c*d^2)/(d*e)^(1/2)*arctan(e*x/(d*e)^(1/2))+1/(a*e^2+c*d^2)*(1/8*e*(a/c)^(1/4)*2^(1/2)*(ln((x^2+(a/
c)^(1/4)*x*2^(1/2)+(a/c)^(1/2))/(x^2-(a/c)^(1/4)*x*2^(1/2)+(a/c)^(1/2)))+2*arctan(2^(1/2)/(a/c)^(1/4)*x+1)+2*a
rctan(2^(1/2)/(a/c)^(1/4)*x-1))+1/8*d/(a/c)^(1/4)*2^(1/2)*(ln((x^2-(a/c)^(1/4)*x*2^(1/2)+(a/c)^(1/2))/(x^2+(a/
c)^(1/4)*x*2^(1/2)+(a/c)^(1/2)))+2*arctan(2^(1/2)/(a/c)^(1/4)*x+1)+2*arctan(2^(1/2)/(a/c)^(1/4)*x-1)))
________________________________________________________________________________________
Maxima [A]
time = 0.50, size = 274, normalized size = 0.81 \begin {gather*} -\frac {\sqrt {d} \arctan \left (\frac {x e^{\frac {1}{2}}}{\sqrt {d}}\right ) e^{\frac {1}{2}}}{c d^{2} + a e^{2}} + \frac {\frac {2 \, \sqrt {2} {\left (\sqrt {a} c d + a \sqrt {c} e\right )} \arctan \left (\frac {\sqrt {2} {\left (2 \, \sqrt {c} x + \sqrt {2} a^{\frac {1}{4}} c^{\frac {1}{4}}\right )}}{2 \, \sqrt {\sqrt {a} \sqrt {c}}}\right )}{\sqrt {a} \sqrt {\sqrt {a} \sqrt {c}} \sqrt {c}} + \frac {2 \, \sqrt {2} {\left (\sqrt {a} c d + a \sqrt {c} e\right )} \arctan \left (\frac {\sqrt {2} {\left (2 \, \sqrt {c} x - \sqrt {2} a^{\frac {1}{4}} c^{\frac {1}{4}}\right )}}{2 \, \sqrt {\sqrt {a} \sqrt {c}}}\right )}{\sqrt {a} \sqrt {\sqrt {a} \sqrt {c}} \sqrt {c}} - \frac {\sqrt {2} {\left (\sqrt {a} c d - a \sqrt {c} e\right )} \log \left (\sqrt {c} x^{2} + \sqrt {2} a^{\frac {1}{4}} c^{\frac {1}{4}} x + \sqrt {a}\right )}{a^{\frac {3}{4}} c^{\frac {3}{4}}} + \frac {\sqrt {2} {\left (\sqrt {a} c d - a \sqrt {c} e\right )} \log \left (\sqrt {c} x^{2} - \sqrt {2} a^{\frac {1}{4}} c^{\frac {1}{4}} x + \sqrt {a}\right )}{a^{\frac {3}{4}} c^{\frac {3}{4}}}}{8 \, {\left (c d^{2} + a e^{2}\right )}} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(x^2/(e*x^2+d)/(c*x^4+a),x, algorithm="maxima")
[Out]
-sqrt(d)*arctan(x*e^(1/2)/sqrt(d))*e^(1/2)/(c*d^2 + a*e^2) + 1/8*(2*sqrt(2)*(sqrt(a)*c*d + a*sqrt(c)*e)*arctan
(1/2*sqrt(2)*(2*sqrt(c)*x + sqrt(2)*a^(1/4)*c^(1/4))/sqrt(sqrt(a)*sqrt(c)))/(sqrt(a)*sqrt(sqrt(a)*sqrt(c))*sqr
t(c)) + 2*sqrt(2)*(sqrt(a)*c*d + a*sqrt(c)*e)*arctan(1/2*sqrt(2)*(2*sqrt(c)*x - sqrt(2)*a^(1/4)*c^(1/4))/sqrt(
sqrt(a)*sqrt(c)))/(sqrt(a)*sqrt(sqrt(a)*sqrt(c))*sqrt(c)) - sqrt(2)*(sqrt(a)*c*d - a*sqrt(c)*e)*log(sqrt(c)*x^
2 + sqrt(2)*a^(1/4)*c^(1/4)*x + sqrt(a))/(a^(3/4)*c^(3/4)) + sqrt(2)*(sqrt(a)*c*d - a*sqrt(c)*e)*log(sqrt(c)*x
^2 - sqrt(2)*a^(1/4)*c^(1/4)*x + sqrt(a))/(a^(3/4)*c^(3/4)))/(c*d^2 + a*e^2)
________________________________________________________________________________________
Fricas [B] Leaf count of result is larger than twice the leaf count of optimal. 1812 vs.
\(2 (245) = 490\).
time = 0.44, size = 3651, normalized size = 10.83 \begin {gather*} \text {Too large to display} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(x^2/(e*x^2+d)/(c*x^4+a),x, algorithm="fricas")
[Out]
[-1/4*((c*d^2 + a*e^2)*sqrt(-(2*d*e + (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2
*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d
^2*e^2 + a^2*e^4))*log(-c*d^2*x + a*x*e^2 + (a*c*d^2*e - a^2*e^3 - (a*c^3*d^5 + 2*a^2*c^2*d^3*e^2 + a^3*c*d*e^
4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^
2*e^6 + a^5*c*e^8)))*sqrt(-(2*d*e + (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e
^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d^2
*e^2 + a^2*e^4))) - (c*d^2 + a*e^2)*sqrt(-(2*d*e + (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*
d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*
d^4 + 2*a*c*d^2*e^2 + a^2*e^4))*log(-c*d^2*x + a*x*e^2 - (a*c*d^2*e - a^2*e^3 - (a*c^3*d^5 + 2*a^2*c^2*d^3*e^2
+ a^3*c*d*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 +
4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))*sqrt(-(2*d*e + (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^
2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*d^
4 + 2*a*c*d^2*e^2 + a^2*e^4))) + (c*d^2 + a*e^2)*sqrt(-(2*d*e - (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2
*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c
*e^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4))*log(-c*d^2*x + a*x*e^2 + (a*c*d^2*e - a^2*e^3 + (a*c^3*d^5 + 2*a^
2*c^2*d^3*e^2 + a^3*c*d*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*
c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))*sqrt(-(2*d*e - (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d
^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e
^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4))) - (c*d^2 + a*e^2)*sqrt(-(2*d*e - (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^
4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^
2*e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4))*log(-c*d^2*x + a*x*e^2 - (a*c*d^2*e - a^2*e^3 + (a*c
^3*d^5 + 2*a^2*c^2*d^3*e^2 + a^3*c*d*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6
*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))*sqrt(-(2*d*e - (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)
*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*
e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4))) - 2*sqrt(-d*e)*log((x^2*e - 2*sqrt(-d*e)*x - d)/(x^2*
e + d)))/(c*d^2 + a*e^2), -1/4*(4*sqrt(d)*arctan(x*e^(1/2)/sqrt(d))*e^(1/2) + (c*d^2 + a*e^2)*sqrt(-(2*d*e + (
c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 +
6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4))*log(-c*d^2*x + a*x*e
^2 + (a*c*d^2*e - a^2*e^3 - (a*c^3*d^5 + 2*a^2*c^2*d^3*e^2 + a^3*c*d*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2
*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))*sqrt(-(2*d*e + (c^
2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*
a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4))) - (c*d^2 + a*e^2)*sqr
t(-(2*d*e + (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c
^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4))*log(-c*
d^2*x + a*x*e^2 - (a*c*d^2*e - a^2*e^3 - (a*c^3*d^5 + 2*a^2*c^2*d^3*e^2 + a^3*c*d*e^4)*sqrt(-(c^2*d^4 - 2*a*c*
d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))*sqrt(
-(2*d*e + (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4
*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4))) + (c*d^2
+ a*e^2)*sqrt(-(2*d*e - (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*
d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*
e^4))*log(-c*d^2*x + a*x*e^2 + (a*c*d^2*e - a^2*e^3 + (a*c^3*d^5 + 2*a^2*c^2*d^3*e^2 + a^3*c*d*e^4)*sqrt(-(c^2
*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c
*e^8)))*sqrt(-(2*d*e - (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2*e^4)/(a*c^5*d^
8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^3*d^4*e^4 + 4*a^4*c^2*d^2*e^6 + a^5*c*e^8)))/(c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^
4))) - (c*d^2 + a*e^2)*sqrt(-(2*d*e - (c^2*d^4 + 2*a*c*d^2*e^2 + a^2*e^4)*sqrt(-(c^2*d^4 - 2*a*c*d^2*e^2 + a^2
*e^4)/(a*c^5*d^8 + 4*a^2*c^4*d^6*e^2 + 6*a^3*c^...
________________________________________________________________________________________
Sympy [F(-1)] Timed out
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Timed out} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(x**2/(e*x**2+d)/(c*x**4+a),x)
[Out]
Timed out
________________________________________________________________________________________
Giac [A]
time = 4.00, size = 336, normalized size = 1.00 \begin {gather*} -\frac {\sqrt {d} \arctan \left (\frac {x e^{\frac {1}{2}}}{\sqrt {d}}\right ) e^{\frac {1}{2}}}{c d^{2} + a e^{2}} + \frac {{\left (\left (a c^{3}\right )^{\frac {1}{4}} a c e + \left (a c^{3}\right )^{\frac {3}{4}} d\right )} \arctan \left (\frac {\sqrt {2} {\left (2 \, x + \sqrt {2} \left (\frac {a}{c}\right )^{\frac {1}{4}}\right )}}{2 \, \left (\frac {a}{c}\right )^{\frac {1}{4}}}\right )}{2 \, {\left (\sqrt {2} a c^{3} d^{2} + \sqrt {2} a^{2} c^{2} e^{2}\right )}} + \frac {{\left (\left (a c^{3}\right )^{\frac {1}{4}} a c e + \left (a c^{3}\right )^{\frac {3}{4}} d\right )} \arctan \left (\frac {\sqrt {2} {\left (2 \, x - \sqrt {2} \left (\frac {a}{c}\right )^{\frac {1}{4}}\right )}}{2 \, \left (\frac {a}{c}\right )^{\frac {1}{4}}}\right )}{2 \, {\left (\sqrt {2} a c^{3} d^{2} + \sqrt {2} a^{2} c^{2} e^{2}\right )}} + \frac {{\left (\left (a c^{3}\right )^{\frac {1}{4}} a c e - \left (a c^{3}\right )^{\frac {3}{4}} d\right )} \log \left (x^{2} + \sqrt {2} x \left (\frac {a}{c}\right )^{\frac {1}{4}} + \sqrt {\frac {a}{c}}\right )}{4 \, {\left (\sqrt {2} a c^{3} d^{2} + \sqrt {2} a^{2} c^{2} e^{2}\right )}} - \frac {{\left (\left (a c^{3}\right )^{\frac {1}{4}} a c e - \left (a c^{3}\right )^{\frac {3}{4}} d\right )} \log \left (x^{2} - \sqrt {2} x \left (\frac {a}{c}\right )^{\frac {1}{4}} + \sqrt {\frac {a}{c}}\right )}{4 \, {\left (\sqrt {2} a c^{3} d^{2} + \sqrt {2} a^{2} c^{2} e^{2}\right )}} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(x^2/(e*x^2+d)/(c*x^4+a),x, algorithm="giac")
[Out]
-sqrt(d)*arctan(x*e^(1/2)/sqrt(d))*e^(1/2)/(c*d^2 + a*e^2) + 1/2*((a*c^3)^(1/4)*a*c*e + (a*c^3)^(3/4)*d)*arcta
n(1/2*sqrt(2)*(2*x + sqrt(2)*(a/c)^(1/4))/(a/c)^(1/4))/(sqrt(2)*a*c^3*d^2 + sqrt(2)*a^2*c^2*e^2) + 1/2*((a*c^3
)^(1/4)*a*c*e + (a*c^3)^(3/4)*d)*arctan(1/2*sqrt(2)*(2*x - sqrt(2)*(a/c)^(1/4))/(a/c)^(1/4))/(sqrt(2)*a*c^3*d^
2 + sqrt(2)*a^2*c^2*e^2) + 1/4*((a*c^3)^(1/4)*a*c*e - (a*c^3)^(3/4)*d)*log(x^2 + sqrt(2)*x*(a/c)^(1/4) + sqrt(
a/c))/(sqrt(2)*a*c^3*d^2 + sqrt(2)*a^2*c^2*e^2) - 1/4*((a*c^3)^(1/4)*a*c*e - (a*c^3)^(3/4)*d)*log(x^2 - sqrt(2
)*x*(a/c)^(1/4) + sqrt(a/c))/(sqrt(2)*a*c^3*d^2 + sqrt(2)*a^2*c^2*e^2)
________________________________________________________________________________________
Mupad [B]
time = 1.59, size = 2500, normalized size = 7.42 \begin {gather*} \text {Too large to display} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(x^2/((a + c*x^4)*(d + e*x^2)),x)
[Out]
(log(a^2*d*e^7 + c^2*d^5*e^3 - c^2*d*x*(-d*e)^(7/2) + 2*a*c*d^3*e^5 + a^2*e^7*x*(-d*e)^(1/2) + 2*a*c*e^3*x*(-d
*e)^(5/2))*(-d*e)^(1/2))/(2*a*e^2 + 2*c*d^2) - atan((((-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/
(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(((-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d
*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(x*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*
a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(512*a^5*c^4*e^9 - 512*a^2*c^7*d^6*e^3 - 512*
a^3*c^6*d^4*e^5 + 512*a^4*c^5*d^2*e^7) + 192*a^4*c^4*d*e^7 + 192*a^2*c^6*d^5*e^3 + 384*a^3*c^5*d^3*e^5) - x*(1
6*a*c^6*d^5*e^2 - 112*a^3*c^4*d*e^6 + 160*a^2*c^5*d^3*e^4))*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c
*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2) + 4*a*c^5*d^4*e^2 + 52*a^2*c^4*d^2*e^4) + x*(2*a
^2*c^3*e^5 - 4*a*c^4*d^2*e^3))*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*
e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*1i - ((-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4
+ a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(((-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*
d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(192*a^4*c^4*d*e^7 - x*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2)
+ 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(512*a^5*c^4*e^9 - 512*a^2*c^7*d^6*e^3 -
512*a^3*c^6*d^4*e^5 + 512*a^4*c^5*d^2*e^7) + 192*a^2*c^6*d^5*e^3 + 384*a^3*c^5*d^3*e^5) + x*(16*a*c^6*d^5*e^2
- 112*a^3*c^4*d*e^6 + 160*a^2*c^5*d^3*e^4))*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^
3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2) + 4*a*c^5*d^4*e^2 + 52*a^2*c^4*d^2*e^4) - x*(2*a^2*c^3*e^5 - 4*
a*c^4*d^2*e^3))*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2
*d^2*e^2)))^(1/2)*1i)/(((-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2
*a^2*c^2*d^2*e^2)))^(1/2)*(((-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4
+ 2*a^2*c^2*d^2*e^2)))^(1/2)*(x*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*
c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(512*a^5*c^4*e^9 - 512*a^2*c^7*d^6*e^3 - 512*a^3*c^6*d^4*e^5 + 512*a^4*c^5*
d^2*e^7) + 192*a^4*c^4*d*e^7 + 192*a^2*c^6*d^5*e^3 + 384*a^3*c^5*d^3*e^5) - x*(16*a*c^6*d^5*e^2 - 112*a^3*c^4*
d*e^6 + 160*a^2*c^5*d^3*e^4))*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e
^4 + 2*a^2*c^2*d^2*e^2)))^(1/2) + 4*a*c^5*d^4*e^2 + 52*a^2*c^4*d^2*e^4) + x*(2*a^2*c^3*e^5 - 4*a*c^4*d^2*e^3))
*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/
2) + ((-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2))
)^(1/2)*(((-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e
^2)))^(1/2)*(192*a^4*c^4*d*e^7 - x*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^
3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(512*a^5*c^4*e^9 - 512*a^2*c^7*d^6*e^3 - 512*a^3*c^6*d^4*e^5 + 512*a^4*c^
5*d^2*e^7) + 192*a^2*c^6*d^5*e^3 + 384*a^3*c^5*d^3*e^5) + x*(16*a*c^6*d^5*e^2 - 112*a^3*c^4*d*e^6 + 160*a^2*c^
5*d^3*e^4))*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2
*e^2)))^(1/2) + 4*a*c^5*d^4*e^2 + 52*a^2*c^4*d^2*e^4) - x*(2*a^2*c^3*e^5 - 4*a*c^4*d^2*e^3))*(-(c*d^2*(-a*c)^(
1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2) + 2*a*c^3*d*e^3
))*(-(c*d^2*(-a*c)^(1/2) - a*e^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(
1/2)*2i - (log(a^2*d*e^7 + c^2*d^5*e^3 + c^2*d*x*(-d*e)^(7/2) + 2*a*c*d^3*e^5 - a^2*e^7*x*(-d*e)^(1/2) - 2*a*c
*e^3*x*(-d*e)^(5/2))*(-d*e)^(1/2))/(2*(a*e^2 + c*d^2)) - atan((((-(a*e^2*(-a*c)^(1/2) - c*d^2*(-a*c)^(1/2) + 2
*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(((-(a*e^2*(-a*c)^(1/2) - c*d^2*(-a*c)^(1/2)
+ 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(x*(-(a*e^2*(-a*c)^(1/2) - c*d^2*(-a*c)^
(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(512*a^5*c^4*e^9 - 512*a^2*c^7*d^6*
e^3 - 512*a^3*c^6*d^4*e^5 + 512*a^4*c^5*d^2*e^7) + 192*a^4*c^4*d*e^7 + 192*a^2*c^6*d^5*e^3 + 384*a^3*c^5*d^3*e
^5) - x*(16*a*c^6*d^5*e^2 - 112*a^3*c^4*d*e^6 + 160*a^2*c^5*d^3*e^4))*(-(a*e^2*(-a*c)^(1/2) - c*d^2*(-a*c)^(1/
2) + 2*a*c*d*e)/(16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2) + 4*a*c^5*d^4*e^2 + 52*a^2*c^4*d^2*e^4
) + x*(2*a^2*c^3*e^5 - 4*a*c^4*d^2*e^3))*(-(a*e^2*(-a*c)^(1/2) - c*d^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(a*c^3*d^
4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*1i - ((-(a*e^2*(-a*c)^(1/2) - c*d^2*(-a*c)^(1/2) + 2*a*c*d*e)/(16*(
a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2)))^(1/2)*(((-(a*e^2*(-a*c)^(1/2) - c*d^2*(-a*c)^(1/2) + 2*a*c*d*e)/(
16*(a*c^3*d^4 + a^3*c*e^4 + 2*a^2*c^2*d^2*e^2))...
________________________________________________________________________________________